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Abstract: We study the effects of CP violation in charged Higgs boson production pp — > 
tH^ 1 + X at the LHC, as well as in the charged Higgs boson decays — > tb and — > 
W^Hf, i = 1,2,3. The study is done in the framework of the type II complex Two 
Higgs Doublet Model (2HDM) with softly broken Zi symmetry. In this model violation of 
CP invariance is induced by the complex parameter mf 2 of the tree-level Higgs potential. 
We calculate the CP violating rate asymmetries for H + and H~ production and decays 
as well as for the combined processes at one-loop level and perform a detailed numerical 
analysis. All calculations are done with the automatic amplitude generator FeynArts and 
the calculational tool FormCalc, for which we have written a complete complex 2HDM 
model file and relevant fortran drivers. The implementation of the complex 2HDM in 
FeynArts and FormCalc is described. In comparison with the analogous results in the 
MSSM, all considered CP violating asymmetries are smaller by an order of magnitude and 
do not exceed 2-^3%. 
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1 Introduction 



The CERN Large Hadron Collider (LHC) has started its operation aiming at a direct 
verification of one of the different candidates that generalize the Standard Model (SM). 
The experiments at the LHC have to distinguish between the predictions of the various 
theoretical models. Having in mind the large number of free parameters that most of these 
models introduce, this task is highly nontrivial. Powerful software and hardware resources 
are required so that their analyses lead to definite predictions for experimental searches. 

An important approach for testing fundamental theories is studying discrete symme- 
tries as model properties. In particular, the CP symmetry is known to be violated in nature 
[1-3]. According to Sakharov's theorem, the mechanism of CP violation (CPV) in the SM 
is not strong enough to explain the baryon asymmetry in the universe - more CPV is 
needed. For that reason all models beyond the SM suggest additional and different sources 
of CPV. 

On the other hand, the Higgs boson is not yet found and the mechanism of electroweak 
symmetry breaking (EWSB) remains the only part of the SM which is not yet verified. The 
extensions of the SM enlarge the Higgs sector and predict both charged and new neutral 
Higgs bosons. If a neutral Higgs boson is discovered at the LHC, there will be a long 
way to determine wether it belongs to the SM or to some of its extensions. A discovery 
of a charged Higgs boson though, would be a clear signal for Physics beyond the SM. A 
possibility to distinguish between the different models containing a charged Higgs is looking 
for effects of CPV in processes with H^, which we shall explore in this paper. 

The simplest extensions of the SM are the models with two Higgs doublets. All of 
them provide new sources of CPV. The most popular one is the Minimal Supersymmetric 
Standard Model (MSSM). The MSSM Higgs sector is a constrained 2HDM of type II. On 
the other hand, the general Complex Two Higgs Doublet Model (C2HDM) has attracted 
much attention due to the CPV it can accommodate [4-6] and due to its simplicity. The 
physical mass eigenstates in the Higgs sector of the 2HDM and the MSSM are the same 
- there are two charged and three neutral Hf, i = 1,2,3 Higgs bosons. Processes 
involving can generate a CPV asymmetry at one-loop level in both MSSM and C2HDM. 
However, the CPV sources in these two models are different. 

In the MSSM, the tree-level Higgs potential is real and thus the neutral Higgs bosons 
have definite CP parities and preserve CP invariance. CPV results from the non-zero CP 
phases of the higgsino mass parameter /x = in the superpotential, the gaugino 

mass parameters Mi = \Mi\ e^', i = 1,2,3, and the trilinear couplings Af = \AA e l ^ f (f 
stands for a fermion) in the soft SUSY breaking Lagrangian. In particular, in the MSSM 
neutral Higgs sector, the presence of these CP violating phases induces mixing between the 
CP-even and the CP-odd scalars at one-loop level, yielding three mass eigenstates [7]. 

In the C2HDM CPV is induced by the complex parameters of the tree-level Higgs 
potential and thus the physical neutral Higgs bosons are mixtures of the CP-even and the 
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CP-odd states. The interactions of these neutral Higgs bosons with fermions and gauge 
bosons violate CP invariance. We stress that in the C2HDM the mixture of the CP-even 
and the CP-odd neutral Higgs bosons is already at tree-level, while in the MSSM it is a 
loop-induced effect, generated by SUSY-loop corrections. 

Effects of CPV in the decays of into "ordinary" (SM) particles in the framework 
of the MSSM were studied in [8—11]. These CPV decay rate asymmetries are of interest 
for a future linear collider, where the charged Higgs will be produced in pairs. There CPV 
occurs only in the decays of H^. For the LHC one must take into account CPV in the 
production process as well. Recently, W^H 11 production at the LHC was studied in [12]. 
In [13-15] we studied CPV in the combined process of a charged Higgs boson production 
at the LHC: 

pp -»• H± t + X, (1.1) 

followed by the subsequent decays of H^, where the production process (1.1) is due to the 
partonic process bg — > t. In both production and decays CPV is induced by one-loop 
radiative corrections with supersymmetric (SUSY) particles in the loops. Our numerical 
analysis showed that the CPV asymmetries, both in the decays and in the production, 
can be rather large. However, in the — > tb decay mode, there can occur cancellations 
reducing the total asymmetry. 

In this paper we perform an analogous study within the complex 2HDM. In addition 
to the production rate asymmetry we calculate the CPV decay rate asymmetries in the 
dominant decay modes of in the C2HDM: 

H^^tb and H ± ^W ± H^, i = l,2. (1.2) 

We also consider the CPV asymmetries in the combined processes of production (1.1) 
and decays (1.2). We work in type II C2HDM with softly broken Z2 symmetry of the 
Lagrangian. In order to perform a numerical analysis, we have generalized the existing 
codes of the FeynArts (FA) and FormCalc (FC) packages [16-19], for calculating processes 
in the C2HDM. 

The paper is organized as follows. In section 2 we shortly review the general 2HDM, 
fix the conventions for the scalar potential, the Yukawa interactions and the parameter set 
we work with. We derive the Higgs couplings to fermions and gauge bosons and list the 
existing theoretical constraints on the 2HDM Lagrangian. In the third section we study 
the production and decay processes mentioned above, present the expressions for the CPV 
asymmetries, and show the one-loop contributions to these asymmetries in the C2HDM. 
We proceed with a detailed numerical analysis in section 4. The numerical results for the 
CPV asymmetries in the charged Higgs decays (1.2), in the charged Higgs production (1.1), 
as well as for the combined process of production (1.1) and subsequent decays (1.2) at LHC 
at the center-of-mass energy y/s = 14 TeV in the C2HDM are presented. We also discuss 
the experimental constraints on the C2HDM parameter space. The numerical analysis 
is done using the FA and FC packages. The implementation of the C2HDM into these 
packages is given in detail in the appendix, which also contains the Lagrangian relevant for 
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our study, the expression for Ap implemented in the code and some important relations 
between the parameters of the scalar potential and the physical Higgs masses. As usual, 
the paper ends with conclusions. 



2 General 2HDM: short review and notations 



2.1 Scalar potential and its parameterization 

The general 2HDM is obtained via extending the SM Higgs sector, consisting of one complex 
Y= +1, SU(2)l doublet scalar field $x> with an additional complex Y = +1, SU(2)l doublet 
scalar field <I>2- Using $1^, one can build the most general renormalizable SU(2)l x U(1)y 
gauge invariant Higgs potential [20-27]: 



V H i ggs (*l,$2) = y(^i: 

l r, 
+ 2 



A2 
2 



+ -^($2$2) 2 + A 3 (*J*i)(*S* 2 ) + A 4 ($ t 1 $ 2 )(^ $ i^ 



A5($i$2) 2 + h.c] + { [A 6 ($t$i) + A 7 ($ 2 $)] (*i*a) + h.c.} 

- - jmfi^f*! + m? 2 $i$2 + h.c. +m| 2 $|$ 2 } . (2.1) 

By hermiticity of eq. (2.1), Ai 2,3,4, as well as my and 777-22 are real- valued; while the 
dimensionless parameters A5, Xq, A7 and mf 2 are in general complex. 

2.1.1 Mass eigenstates 

After the SU(2)l x U(1)y gauge symmetry is broken down to U(l) em via the Higgs mecha- 
nism, one can choose a basis where the vacuum expectation values (VEVs) of the two Higgs 
doublets, v\ and t> 2 are non-zero, real and positive, and fix the following parameterization 
[28, 29]: 



$1 



(j, _ [ ft \ (2 21 

(vi + m + iXi)/V2 J ' 2 \(v 2 + V2 + iX2)/V2 J ' k ' ' 

Here 771,2 and xi,2 are neutral scalar fields and y?^ 2 are charged scalar fields. The physical 
Higgs eigenstates are obtained as follows. 

The charged Higgs fields and the charged would-be Goldstone boson fields are 
a mixture of the charged components of the Higgs doublets (2.2), (ff 2 : 



— sin fiipf + cos (3(pf , 

cos f3(pf + sin (3tpf , (2-3) 

where the mixing angle /3 is defined through the ratio of the VEVs of the two Higgs doublets 
$2 and 3>i, tan/3 = u 2 /ui. give masses to the bosons. 

Obtaining the neutral physical Higgs states is a few steps procedure. First, one rotates 
the imaginary parts of the neutral components of eq. (2.2): (xi> X2) into the basis (G°, 773): 1 

G° = cos/?xi + sin/3x2 , 

773 = - sin /?xi + cos /?X2 , (2.4) 



x Note that in the case of nii 2 = A6 = A7 = and all other parameters of eq. (2.1) are real, the physical 
Higgs sector of the 2HDM is analogous to the one of the tree-level MSSM. In this case the scalar field 773 is 
equivalent to the MSSM neutral CP-odd Higgs boson A . 
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where G° is the would-be Goldstone boson which gives a mass to the Z gauge boson. After 
elimination of the Goldstone mode, the remaining neutral CP-odd component 773 mixes 
with the neutral CP-even components 77^2 ■ The relevant squared mass matrix -Mfj = 
9 2 ^Higgs/ (drjidrjj), i,j = 1, 2, 3, has to be rotated from the so called "weak basis" (771, 772, 773) 
to the diagonal basis (Hf, H®, H®) by an orthogonal 3x3 matrix TZ as follows: 

KM 2 K T = .Mdiag = diag(M^ ? , M^o, M^o) , (2.5) 

with 

(2.6) 

where we have defined the Higgs fields Hf such that their masses satisfy the inequalities: 

M H o < M H o < M H o . (2.7) 

Note that the mass eigenstates Hf have a mixed CP structure. 

Following [30] , we parametrize the orthogonal 3x3 matrix TZ by three rotation angles 
a,:, i = 1, 2, 3: 

10 \ / cos 0:2 sina^A / cosai sinai 
TZ = I cosa3 sina3 I I 1 I I — sinai cosai 
— sin a3 cos «3 / V — sin a 2 cos a% I \ 1 

-(ci s 2 s 3 + si c 3 ) C1C3-S1S2S3 c 2 s 3 , (2.8) 
-C1S2C3+S1S3 -(ci S3 + si s 2 c 3 ) c 2 c 3 / 

with Sj = sinaj and a = cosaj, which we vary in our numerical analysis in the following 
ranges: 

TT IT TT TT TT , 

"2<«i<2 ; ~2 <a2 ~2 ; °^ a3 ^2' (2 ' 9) 
2.2 Z2 symmetry and input parameter set 

In the most general 2HDM, some types of Yukawa interactions can introduce flavour chang- 
ing neutral currents (FCNC) already at tree level. It is well known that the latter effects 
are small in nature. This problem has been solved by imposing a discrete Z2 symmetry 
on the Lagrangian. It forbids $1 -B- $2 transitions and in its exact form it also leads to 
conservation of CP [31]. In order to allow some effects of CPV it is necessary to violate the 
Z2 symmetry. Basically, there are two ways of Z2 symmetry violation - "soft" and "hard". 
A softly broken Z2 symmetry suppresses FCNC at tree level, but still allows CPV. 

In this paper we will work in a model of a softly broken Z2 symmetry of the 2HDM 
Lagrangian. This forbids the quartic terms proportional to Xq and A7 in eq. (2.1), but the 
quadratic term with m\ 2 is still allowed [32, 33]: 

= y (<^i) 2 + y(*!*2) 2 + A 3 ($ t 1 <i>i)(<I> 2 <I> 2 ) + A^I^X^i) 



A 5 ($ t 1 $ 2 ) 2 + h.c. - ^|m 2 1 $|^>i+ m 2 2 $l$ 2 + h.c. +m^$^ 2 }- (2.10) 
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The Higgs potential (2.10) has 12 real parameters: 2 real masses: 22 , 2 VEVs: 
Ul,2) four real quartic couplings: Ai )2 ,3,4 and two complex parameters: A5 and m 2 2 . The 
conditions for having an extremum of eq. (2.10) reduce the number of parameters: m 2 ^ 22 
are eliminated by the minimization conditions, and the combination v\ + v\ is fixed at the 
electroweak scale v = (y^G^) -1 / 2 = 246 GeV. Moreover, in this case the minimization 
conditions also relates Im(m| 2 ) and Im(As): 

Im (m 2 2 ) = f 1 f 2 Im (A5) . (2-11) 

Thus, our Higgs potential (2.10) is a function of 8 real independent parameters: 

{Ai,2, 3 ,4, Re (A 5 ), Re(m 2 2 ), tan/3, Im(m 2 2 )}. (2.12) 

It contains minimal CPV generated by m 2 2 7^ and complex. In our further analysis we 
will use the following parameter set equivalent to eq. (2.12): 

|m#o, M h o, M h +, ax, a 2 , 0-3, tan/3, Re(mi 2 )|. (2-13) 

Note that the mass of the heaviest neutral Higgs boson H® is not an independent parameter. 
In the considered CP violating case, the matrix elements (.M 2 )i3 and (.M 2 )23 of the squared 
mass matrix (2.5) are non-zero and correlated [34]: 

(.M 2 ) 13 = tan/3 (_M 2 ) 23 . (2.14) 
Writing this relation in terms of the physical masses M^o, M^o, M^o one obtains [34]: 

2 M 2 H0 R 13 (R 12 t&nf3 -Rn) + M^ R 23 (R22 tan /3 - R 2l ) 
Mh °* = ~ i233(i23i-i232tan/3) ' (2 " 15) 

where Rij ,i,j = 1, 2, 3, are the elements of the rotation matrix (2.8). 

The expressions for the parameters Ai.,2,3,4, ReA5,ImA5 of the scalar potential (2.10) 
as functions of the physical masses and mixing angles are given in appendix B. 

In appendix C we also list the triple scalar couplings of the Higgs bosons calculated from 
the potential (2.10) and relevant to our study. Note, that the couplings C(H° H ± G T ), i = 
1, 2, 3, have both real and imaginary parts and can lead to CPV. 

2.3 Higgs and gauge boson interactions 

The Higgs and gauge boson interactions arise from the covariant derivatives of the doublet 
fields <I>i 2: 



-gauge — 

i=l 



^2(D^(D^i), D" = d" + igf a W£ + ig'YiB^/2 (2.16) 



where T a are the isospin generators, Yi are the hypercharges of the Higgs and 
are the SU(2)l and U(l)y gauge fields, g, g' are the corresponding SU(2)l and U(l)y gauge 
couplings. 
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After having rotated the Higgs and gauge bosons fields to their mass eigenstates bases 
one obtains terms of triple and quartic interactions between them. The interactions relevant 
to our study are HjWW and H? tP, with the corresponding couplings 

C(HjWW) = cos PRa + sin fiR i2 , 
C(H° W ± H T ) = T i{smPRn - cos (3R i2 ) ± R i3 . (2.17) 

Note that C(H?WW) is purely real, while C(H? IP) has both real and imaginary parts 
and can lead to effects of CPV, as we will see later. The relevant Lagrangian is given in 
appendix C. 

From the condition for unitarity of the rotation matrix 1Z one derives the following 
sum rules: 

C(H°WW) 2 + |C(H° W + IT)| 2 = 1 for each i = 1,2,3 (2.18) 

3 

J^C(H°WW) 2 = 1 (2.19) 

i=i 

From eqs. (2.18) and (2.17) follows that if for a fixed i the term |C(H£ W+ H~)| 2 is sup- 
pressed, then (sin/3i?ji — cos (3Ri 2 ) 2 ~ and i? 2 3 « 0. The relation Rj 3 « implies that il? 
has a very small pseudoscalar component and is dominantly a CP-even state. Furthermore, 
the sum rule (2.19) implies that the other Higgs states decouple, i.e. C(H°WW) 2 « for 
j ^ i- We will come back to the physical consequences of these sum rules in our physics 
discussion on the studied decay and production processes in section 4. 

2.4 Yukawa interactions 

In the framework of the 2HDM various models of Yukawa interactions can be realized [4]. 
Depending on the Yukawa interaction one distinguishes between different types of 2HDM's. 

The most general 2HDM, in which each fermion doublet and singlet couples to both 
Higgs doublets is called the type-Ill model. This model leads to FCNC already at tree 
level and is a subject of severe constraints from flavour physics observables [21, 35]. As we 
already discussed, in order to avoid problems with FCNC usually a Z 2 discrete symmetry 
is imposed on the Lagrangian [31]. 

In the 2HDM type-I all fermions, quarks and leptons, couple to only one of the Higgs 
doublets exactly like in the SM. These models are interesting as the decoupling Higgs is a 
natural candidate for dark matter. 

There exists another type of 2HDM where one of the Higgs doublets couples to leptons 
and the other one to the up and down quarks [36-38]. The phenomenology of these models 
has been reviewed recently in [39-41]. 

In the present paper we work in type-II 2HDM. In this model, down-type quarks and 
charged leptons couple to <&i and the up- type quarks couple to the other Higgs doublet <3?2- 
( In the MSSM the Higgs sector is also a 2HDM type-II, but with the Higgs self-interactions 
fixed by gauge couplings. ) In order to avoid FCNC at tree level, but allow for CPV, our 
discrete Z 2 symmetry is softly broken, i.e. A6 = A7 = 0, but m\ 2 is non-zero and complex, 
see section 2.2. 
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The parts of the interaction Lagrangian of Higgs bosons and fermions relevant for our 
discussions are given in appendix C. Since the neutral Higgses are mixtures of CP-odd and 
CP-even states, their couplings to a fermion pair have the general form a + % b 75 with a 
and b real, that can lead to CPV. 



2.5 Theoretical constraints on the Higgs potential 

The Higgs potential given by eq. (2.10) has to satisfy some general requirements 

like positivity, unitarity and perturbativity [30]. These requirements together with the 
minimum conditions naturally lead to constraints on its parameters [28, 30, 34]. The 
theoretical constraints on the potential (2.10) are well studied and will not be a subject 
of this note. For completeness, we list below the expressions we have implemented in our 
numerical code. 

In order to have a stable vacuum the potential should be positive for large values of 
\4>k\, which leads to the constraints [28]: 



Ai>0, A 2 >0, A 3 + VA1A2 > 0, A 4 + A 4 -|A 5 | + v / ATA 2 >0. 



(2.20) 



From the requirement that theory remains perturbative we have: |Aj| < 8n for any i. The 
unitarity requirement means that the tree-level amplitudes for Higgs-Higgs, Higgs-vector 
boson and vector boson- vector boson scattering should not exceed the unitarity limit when 
contributing to the s-wave. Both unitarity and perturbativity requirements lead to the 
constraints [28]: 



1 



Ai + A 2 ± V(Ai-A 2 ) 2 + 4|A 5 | : 



Ai + A 2 ± 



+ 4AJ 



3(Ai + A 2 ) ± V / 9(A 1 -A^ + 4(A 3 + A 4 ) 2 ) 



A 3 ±A 4 |<87r, |A 3 ±|A 5 



2 

< 8?r, 



< 8vr, 



< 8vr, 



< 8vr, 



A3 +2A 4 ± I A 5 1 1 < 87T. 



(2.21) 



As the eqs. (2.20) and (2.21) only depend on the absolute value of A5, they do not constrain 
the phase of A5. We will see later, that these theoretical constraints already strongly reduce 
the C2HDM parameter space [28, 30, 34]. 



3 CP violation in H ± production and decays 
3.1 The processes 

We study CPV induced by one-loop corrections in the C2HDM in the following processes 
involving the charged Higgs boson: 
• Decays: 

H^-ttb, H ± ^W ± H i , * = 1,2. (3.1) 
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• Associated production with a top quark at the LHC: 

pp -> HH + X , (3.2) 

with the partonic subprocesses: 

b g^t H~ and b g ^ t H + . (3.3) 

• Charged Higgs production (3.2) plus subsequent decays (3.1). 

The tree- level graphs of the considered processes for H~ are shown in figure 1. 




Figure 1. The tree- level graphs for bg — > tH : a) s-channel, b) t-channel, and charged Higgs 
decays: c) H~ ib, d) H~ W~H?. 



3.2 CP violating asymmetries 

For the processes listed in section (3.1) we investigate the following CP violating rate 
asymmetries: 

• Decay rate asymmetries A*^, defined by: 

A DJ i H -»•/)- 2r tree^+ ^ ft ' ( 3A ) 

where D stands for decay and / stands for any of the decay modes: / = tb; W H® with 
i = 1,2. 

• Production rate asymmetry Ap P , related only to the process (3.2) and defined by: 

CP _ a{pp -» H+t) - a(pp -» H~t) 
Ap ~ 2a^{pp^H+t) ' [6 - b) 

where P stands for production. 

• Asymmetries A^ p for the combined processes of production and subsequent decays, 
defined by: 



a cp = °(PP -> tH+ -> if) -_v(pp -+tH -» tf) i ,, 



In [13] we have shown that in narrow width approximation, when the decay width of H + 
is much smaller than its mass, the total asymmetry A^j p given by eq. (3.6) is an algebraic 
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sum of the asymmetry A^, p in the production and the asymmetry A^ P j in the decay / of 
the charged Higgs boson: 



A c f p = A C P P + A% p . 



(3.7) 



The decay asymmetries (3.4) might be of interest for the ILC, where CPV can occur 
only in the decays. Since the measurements on b — > 57 put a stringent lower limit on 
the charged Higgs mass M H ± > 295 GeV, the decay modes (3.1) are dominant. 

3.3 CP violating loop contributions 

In order to get CPV using the asymmetries introduced in section 3.2, we need both non- 
zero CP violating phases in the Lagrangian and CP conserving phases (strong phases) in 
the absorptive parts of the one-loop amplitudes. 



H 




H 




11 






Figure 2. Generic CP violating selfenergy and vertex contributions to H — ► tb. F denotes a 
generic fermion, S a generic neutral or charged scalar field, and V a generic vector boson. 



In the C2HDM, the CP violating phases arise from: 

• Neutral Higgs couplings to a fermion pair 

• Charged Higgs - Neutral Higgs - gauge bosons couplings 

• Charged Higgs - Neutral Higgs - Goldstone bosons couplings 

The parts of the C2HDM Lagragian needed for this study are given in the appendix C. 

The CP conserving phases originate from various on-shell intermediate states of the 
one-loop amplitudes. For the considered production and decay processes, the strong 
phases coming from cuts e.g. in t — > bW^ and H® — > bb will always contribute, while 
the strong phases coming from cuts in — > W^Hf, — > G^Hf, — > tb, etc. in 
the loops, will contribute only if they are kinematically allowed. All possible generic CP 



- 10 - 



„0 ttO u tj-0 ttO 

til til til til til 




U V 



Figure 3. Generic CP violating selfenergy and vertex contributions to H —> W H®. 



violating one-loop contributions to the decay — > tb are shown in figure 2 and to the 
decay -»■ W ± Hf in fi gure ,3. 

Furthermore, we have three types of possible generic CP violating loop contributions 
to the partonic cross sections of the production processes (1.1): selfenergy contributions, 
shown in figure 4; vertex contributions shown in figure 5 and box diagram contributions 
shown in figure 6. All diagrams have been generated with FA [16] package, for which we 
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have written a complete model file for the C2HDM, see Appendices D.l and D.2 for details. 

f 




f g h 



Figure 4. Generic CP violating selfenergy contributions to charged Higgs boson production b g — >• 
H+l 



t t t t b F 




a b c d e 



b $ b S b V b b 




f q h k I 



Figure 5. Generic CP violating vertex contributions to charged Higgs boson production b g — > H + t. 




a b c d e 



Figure 6. Generic CP violating box contributions to charged Higgs boson production b g — » H + t. 
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4 Numerical results 



In this section we present our numerical results for the CP violating asymmetries (3.4), 
(3.5) and (3.6) in the C2HDM with a softly broken Z2 symmetry (2.10). All calculations 
have been done using the packages FA and FC [16], for which we have written a complete 
FA model file and have extended the corresponding FC fortran drivers for the C2HDM. 
The implementation is described for the FA model file in the appendix D.l and for the 
FC fortran drivers in D.2. In our numerical analysis we also have used LoopTools [42-44]. 
The calculations are done in the 'tHooft-Feynman gauge using dimensional regularization. 
However, we stress that we do not renormalize any parameters or fields since our CPV rate 
asymmetries involve only the imaginary parts of the loop integrals which are always finite. 
For the evaluation of the PDF's we use CTEQ5L, with a s , calculated at the scale Q = 

In this paper we work with the following set of real input parameters [30] , see section 

2.2: 

|m^o, M h o, M h +, a>i, ai2, «3, tan/3, Re(mi 2 )l. (4.1) 

In the literature often the parameter fi is used instead of Re (77112): 

v 2 

f j 2 = - Re(m? 2 ), v 2 = v\ + v\ . (4.2) 

2tjit;2 

The expressions for the parameters of the scalar potential (2.10), Ai, i = 1,2,3,4,5, as 
functions of the physical masses, mixing angles and \i [34] are given in appendix B. For the 
values of the SM parameters used we refer to [45], except for the top mass, taken from the 
last Tevatron measurement [46]: 

mt = 173.1 GeV , m b = 4.7GeV, and a = 1/137.03599. (4.3) 

Furthermore, we would like to add some comments on the existing constraints on the 
values of the C2HDM parameters. Basically they come from: 

• theory - these are the requirements for positivity and unitarity of the Higgs potential, 
see section 2.5 

• experiment - mainly coming from the electroweak precision data at LEP. These 
constraints we will discuss in the next subsection 

4.1 Experimental constraints 

In principle, there is quite a long list of experimental bounds to constrain the phenomenol- 
ogy of the 2HDM. Examples are: B-mixing/B-decays constraints, LEP2 non-discovery 
constraints, direct searches for H , muon anomalous magnetic moment, electron electric 
dipole moment (EDM), etc., see [30, 33, 47, 48] and the references therein. Many of these 
constraints are strongly dependent on the model of Yukawa interactions, on the consid- 
ered process and on the related parameter space. As already mentioned, the theoretical 
constraints described in section 2.5 already strongly reduce the parameter space of the 
C2HDM [48] . In [35] the authors combine these constraints with the existing experimental 
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constraints and review the "profile of the surviving parameter space" of the type II 2HDM. 
There they show that large values of tan/3 are forbidden by the unitarity constraints, see 
section 2.5, except for the case M H o < p [35]. Furthermore, they show that the experi- 
mental constraints also usually exclude parameter regions for tan/3 ~ 10 [35]. However, 
in principle it can be shown that one can fine-tune the parameters to allow some tiny 
parameter regions for large tan/3 [49]. 

Among the experimental constraints we consider the following ones, which have the 
strongest impact on restricting the general 2HDM parameter space: 

• For the lightest neutral Higgs boson we take into account the LEP2 non-discovery 
bound, M H o > 114.4 GeV [50]. 

• In addition to the lower bounds from the CERN LEP and Tevatron direct searches 
[51-54], the charged Higgs mass is quite severely constrained by the B — > X s ^ data [55-61]. 
At the next-to-next-to-leading order in the type-II 2HDM it implies that M H ± > 295 GeV, 
[59-61]. In our analysis we follow the latter lower bound. 

• Another important experimental constraint, coming from the electroweak physics, is 
related to the precise determination of the p-parameter [62]. Its deviation Ap from the SM 
value should accommodate all new physics contributions. Ap is constrained by the error 
of the measured value of the parameter po, which at the 2 a level is [50]: 

po = 1.0004 t°oS ■ (4-4) 

In our numerical code we implement the expressions for Ap given in [30], with the require- 
ment: 

-0.0011 < Ap < 0.0029. (4.5) 

The analytic expressions for the extra contributions in Ap in the framework of the C2HDM 
are given in appendix A. 

4.2 CP violation in i7 ± -decays 

We discuss the CPV decay rate asymmetries Ap P p given with eq. (3.4), for the decays 
E ± tb and -> W ± H?, i = 1,2. Note that we do not consider the decay mode 

— > W^H® as a possible measurable channel for the CPV asymmetry due to its tiny 
branching ratio (BR) in the considered parameter ranges. 

In figure 7, figure 8 and figure 9 we present the CP asymmetries , A^ wh q 

and A^ p wh0 , respectively, as functions of the charged Higgs mass, for several values of 
tan/3. The curves in the figures are cut once any of the parameter constraints discussed in 
section 2.5 and section 4.1 is violated. In particular, in the considered parameter range, the 
Ap constraint imposes that the charged Higgs mass Mjj+ cannot acquire large values: in 
figure 7, for tan /3 = 1.5, M H + cannot exceed the limit of ~ 368 GeV, while for tan /3 = 4.5: 
M H+ < 408 GeV, etc.. 

All discussed CPV asymmetries exhibit a mild dependence on M H + and a strong 
dependence on tan/3: Ap P tb increases, while A^ wu(i and A^^^ strongly decrease when 
tan/3 increases. This is shown in figures 7, 8, and 9. 
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Figure 7. Left: CPV asymmetry Ap P tb as a function of M H + for four values of tan/3. The other 
parameters are: M H o = 120 GeV, M H o = 220 GeV, Re(mi 2 ) = 170 GeV, a x = 0.8, a 2 = -0.9 and 
a 3 = 7r/3. Right: CPV vertex and selfenergy contributions to Af) P tb , as functions of the charged 
Higgs mass for tan/3 = 1.5. 
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Figure 8. Left: CPV asymmetry A ( ^ ) P WH0 as a function of M H + for four values of tan/3. The 
other parameters are the same as for figure 7. Right: CPV vertex and selfenergy contributions to 
■^dwh^ as f unc ti° ns of the charged Higgs mass for tan/3 = 1.5. 

As we already commented in section 4.1 large values of tan/3 are excluded, except for 
some small areas in the region M H o < [i which require fine-tuning of the parameters [35]. 
Moreover, we have checked that very often for large tan (5 the mass of the heaviest Higgs 
boson M H o becomes tachyonic. Therefore in most of the cases we show numerical results 
for tan/3 in the range 1.5 -r 4.5. 

In figure 7 we see that the asymmetry A^fL is positive and can reach ~ 2.5% for 
tan/3 = 4.5 and M H + ~ 400 GeV. In the same figure we show the vertex and the selfenergy 
contributions in Ap P tb for tan /3 = 1.5. The figure shows that these contributions enter with 
opposite signs and there is a partial cancellation between them. In contrast to the MSSM, 
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Figure 9. CPV asymmetry A c ^ ) P WHa (3.4), as a function of M H + for four values of tan j3. The other 
parameters are: M H o = 120 GeV, M H o = 220 GeV, Rc(to 12 ) = 170 GeV, a x = 0.8, a 2 = -0.9 and 
a 3 = 7r/3. 



here the dominant contribution in the CPV asymmetry comes from the vertex diagrams. 

For the bosonic channel — > W^H®, the CPV asymmetry is shown in figure 8. It is 
seen that ^jf WH o is negative and can reach ~ 1.7% for tan/3 = 1.5 and M H + close to the 
lower limit ~ 297 GeV. In this channel there is also a cancellation between the selfenergy 
and the vertex contributions, which is shown in the figure. In the CPV decay asymmetry 
^f WH o the dominant contribution stems from the selfenergies. 

In figure 9 we show the asymmetry A^^^ as a function of the charged Higgs mass. 
In contrast to the other bosonic channel — > WHf, it is positive and the absolute value 
is larger. For tan (3 = 1.5, it is almost constant ~3.7%. For this asymmetry we do not show 
explicitly the individual contributions. However, we would like to note, that cancellation 
occurs between the vertex and the selfenergy contributions. Similar to the decay into 
the selfenergy contributions are dominant in A ( ^' wh0 . 

In all three discussed cases we have checked that the influence of the CKM matrix is 
very small, therefore we work with a diagonal one. 

In figure 10 we show a scan over the parameters a\ and 02, m order to illustrate 
the allowed domain for (01,02) together with the size of the CPV asymmetry in the 

— > tb decay. As we already mentioned, the excluded parameter regions are due to the 
combination of all theoretical and experimental constraints: vacuum stability, unitarity, 
large additional contributions to Ap, as well as tachyonic modes of M H o or violation of 
the ordering M H o < M H o < M H o . In the figure we see that for smaller values of tan (3 the 
allowed domain is larger and almost vanishes for tan/3 = 4.5 and that the CPV asymmetry 
■^Dtb excee ds 2% only in some small areas. 

In figure 11 we show a similar scan for the — > W^H® decay. When we scan 
over ct\ and 02 for a fixed value of 0:3, we find a region where the coupling C^^H^W^) 
becomes very small. It is seen in the figure that for tan/3 = 1.5 C(H^H^W^) — > for 
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Figure 10. The allowed parameter regions in the (0:1,0:2) plan in the C2HDM together with the 
absolute value of the CPV asymmetry A% p tb . We have taken M H o = 120 GeV, M H o = 220 GeV, 
M H ± = 350 GeV, Re(m 12 ) = 170 GeV, and a 3 = tt/3. On the top left plot tan/3 = 1.5, top right: 
tan/3 = 2, down left: tan/3 = 3, and down right: tan/3 = 4.5. 

«2 ~ and ai ~ tt/3 (yellow color). The position of this region shifts to the right for 
larger values of tan/3. The coupling C(HiH^W^) is complex. Therefore, the condition 
CiH^H^WT) -)• also implies that its imaginary part R13 — > 0, which furthermore means 
that H® is dominated by its CP-even component. 

This can be seen also from the sum rule given by eq. (2.18), which tell us that if 
\C(H^H + W-)\ 2 -)• then C(H^WW) 2 -)• 1 and which has a maximal coupling to a 
W pair, is dominated by its CP-even component. According to our definitions of the CPV 
asymmetries in section 3.2, when ^(H^H^W^)] 2 — > 0, the tree-level width p trec (ff± _). 
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Figure 11. The allowed parameter regions in the (ai,a2) plan in the C2HDM together with the 
absolute value of the CPV asymmetry Aq WhQ . The other parameters are the same as for figure 
10. 



W^Hf) — > and consequently the CPV asymmetry will increase considerably. However, 
the large rate asymmetry we would obtain in this case, would go together with a small BR 
of the i2"± W ± Hl decay. 

In figure 12 we show the BR of H + — > W + H® as a function of a\ for various values of 
tan/3 and the other parameters fixed as in figure 10. It can be seen in the plots that in some 
cases the BR of H + — > W + H® can be larger than 80% together with a CPV asymmetry 
of a few percent. 
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Figure 12. The BR(H + — > VK + i/°) as a function of ai with a2 in the allowed parameter range. 
The other parameters are the same as for figure 10. 

4.3 CP violation in //^-production 

We also study the CPV asymmetry A^ p , given by eq. (3.5), for charged Higgs production 
(3.2) at the LHC with y/s = 14 TeV. For a direct comparison with the results of section 
4.2, in figure 13 we present the asymmetry Ap as a function of the charged Higgs mass 
M H + for the same parameter set as used for figure 7. One can see that here the asymmetry 
is of the same order of magnitude or smaller than in the — y tb. The absolute value 
increases with tan f3, and can go up to 2% for tan /3 = 4 and M H + = 400 GeV. But in 
contrast to the — > tb decay the production asymmetry is negative. Therefore, the total 
asymmetry Ap^ p + Ap P b is very small. However, Ap^ p has the same sign like A^ h0 and 
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according to eq. (3.7) this will increase the total asymmetry in these bosonic modes. 
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Figure 13. CPV asymmetry A c p p (3.5); Left: function of charged Higgs mass for four values of 
tan j3 and the other parameters are as in figure 7: M H a — 120 GeV, M H o = 220 GeV, Re(mi2) = 170 
GeV, a\ = 0.8, a-i — —0.9 and — 7r/3; Right: vertex, selfenergy and box contributions to Ap P 
as a function of Mjj+ for tan j3 = 1.5. 

In figure 13 we show the vertex, selfenergy and box contributions to A^, p for tan/3 = 

I. 5. The vertex contribution is negative and dominating, especially for smaller masses of 
the charged Higgs. The selfenergies and box contributions are both positive, but their sum 
partially cancels with the contributions of the vertex diagrams. 

In figure 14 we show a scan over the angles a\ and a2, together with the absolute value 
of the asymmetry Ap P , as for the decays in the previous section, see figure 10 and figure 

II. The plots are very similar in size and the asymmetry can reach ~ 2%. 

In order to have a better overview on the CPV asymmetry in a larger parameter 
space and having in mind that the integration over the PDFs slows down the calculation 
considerably, we perform a numerical scan over the C2HDM parameters (2.13) using Grid 
computing: 

M H o = 115-f- 125 GeV, with step size 5 GeV , 
M H o = 150 -=- 400 GeV, with step size 50 GeV , 
M H + = 300 -=- 550 GeV, with step size 25 GeV , 
Re (77112) = 10 4- 460, with step size 50 , 
tan (3 = 1 -7- 8, with step size 1 , 

ot,\ = tt/2 -t- vr/2, with step size tt/9 , 

CK2 = tt/2 -t- vr/2, with step size 7r/9 , 

03 = 4- 7r/2, with step size 7r/9 , 

(4.6) 
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Figure 14. The allowed parameter regions in the (0:1,0:2) plan in the C2HDM together with the 
absolute value of the CPV asymmetry Ap P . The other parameters are the same as in figure 10. 



The total number of the scanned parameter points is quite large: 6415200. After 
having applied the theoretical and experimental constraints, see section 2.5 and section 
4.1, the allowed number of parameter points relegates to 67861, which is only 1% of the 
scanned parameter space. 

The numerical results based on the parameter scan (4.6) show that in most of the cases 
the asymmetry is very small, practically zero. The non-zero asymmetry is distributed in 
two non-symmetric bunches - to the right (up to ~ 2%) and to the left (down to ~ —3%) 
from the zero. The non-symmetric structure of the distribution is due to the asymmetric 
Ap constraints, see section 4.1. Furthermore, the results from the numerical scan (4.6) 
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Figure 15. The maximal value of the asymmetry |Ap P | as a function of ot\ and oli - left; and 
M H + and tan/3 - right, based on the parameter scan (4.6). 

show that the CPV asymmetry Ap P can reach 3.5 % (negative) for single isolated points 
of the parameter space. 

In figure 15 we show the maximal value of |^4p P | as a function of a% and «2 - left; 
and M H + and tan/3 - right, based on the parameter scan (4.6). One sees small regions 
where the asymmetry can go up to ~ 3%. These are roughly: {a\,— 0:2) £ (tt/6 + 7r/3) 
and relatively large Mg+ and tan/3: M#+ E (460 -j- 550) and tan/3 E (2-j-4). However, 
our experience showed that making a plot in such region is a question of fine-tuning of the 
parameters due to the severe theoretical constraints ( see section 2.5 ), which often cut the 
parameter space into unconnected subspaces. 

4.4 CP violation in //^-production and subsequent decays 

Combining the production process (3.3) with the subsequent decay — > tb, in [13] 
we have shown analytically that the charged Higgs selfenergy contributions to the total 
asymmetry from the production and the decay exactly cancel. Moreover, we have shown 
that in the MSSM the vertex contributions from the production and from the decay also 
partially cancel. Eventually, the main contribution in the total asymmetry comes from the 
MSSM box contributions to the production. 

In this paper we study the same production process, but in the C2HDM. The selfenergy 
cancellations observed in [13] are not model dependent and we expect them to occur again. 
In figure 16 we show the asymmetry Af b p as a function of the charged Higgs mass. One 
sees that the total asymmetry is of an order of magnitude smaller in comparison to the 
individual decay and production asymmetries. The cancellations are easily traced back 
comparing the selfenergy and vertex contributions in figure 7 and figure 13, which are of 
the same magnitude, but with opposite signs. However, no such cancellations occur in the 
bosonic modes: — > W^H^ 2 - I n figure 17 we show the asymmetries Ayf H and Ay^ H2 



-22- 



T 

Q. 
Q. 

CL 

O 
< 




0.002 



295 310 



330 350 370 
m H+ (GeV) 



390 410 



Figure 16. CPV asymmetry A^- p (3.6) as a function of the charged Higgs mass for four values of 
tan j3 and the other parameters are as in figure 7: M H o = 120 GeV, M H o = 220 GeV, Re(mi2) = 170 
GeV, a\ = 0.8, ot2 = —0.9 and — n/3. 
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Figure 17. CPV asymmetries Ayf Hi and A^f H2 (3.6) as functions of M H + for four values of 
tan/3. The other parameters are: m H a — 120 GeV, m H o — 220 GeV, mi 2 = 170 GeV, ai = 0.8, 
oi2 = —0.9, and — 7r/3. 

(3.6) as functions of Mjj+ for four values of tan /3. As one can see, here the effect is bigger 
and the asymmetry Ay^ Hi can reach up to 2.5% (negative), while the asymmetry Ayf H2 
can go up to 3% (positive). 

In the analogous study performed in the MSSM [13] the main contribution for the 
case with the decay — > tb is due to box graphs with gluino exchange which involves 
the strong coupling constant. This asymmetry reaches its largest value of ~ -12% for 
m H + ~ 550 GeV with a branching ratio BR(H + ->• tb) ~ 20%. 

In contrast to [13], in the present study the asymmetries are due to the exchange of 
neutral Higgs bosons and therefore of electroweak nature. In the tb mode the asymmetry 
always remains below ~ 1% with BR(#+ -)• tb) ~ 50% for m H + ~ 350 GeV. In the W ± H l 
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mode it can reach 2.5% for m H + ~ 350 GeV with BR{H + -> W ± H l ) up to ~ 45%. 



The production rate of H + at the LHC for = 350 GeV and tan j3 = 1.5 is ~ 450 fb 
(including a K- factor from QCD of 1.3, see [64]). Assuming a branching ratio of 50% and 
an integrated luminosity of 200 fb -1 we get N ~ 45000 and y/N ~ 212. Furthermore, 
assuming A CP ~ 0.02, the statistical significance \fNA^ b p ~ 4.2. However, in a realistic 
study the actual signal production rate will be most likely reduced. Moreover, the total 
background rate for charged Higgs production at the LHC is quite large [64]. The obtained 
statistical significance might be too low for a clear observation in the first stage of the LHC 
and therefore an upgrade of its luminosity would be necessary to await. 

5 Conclusions 

In the type II complex 2HDM with softly broken Z2 symmetry, the non-zero and complex 
mf 2 parameter of the tree-level Higgs potential gives rise to CP violation in the production 
process pp — > H^t + X, and in the dominant decay modes of to tb, and to WHi, 
i=l,2. We have calculated the corresponding CP violating rate asymmetries at one- loop 
level both in the production and in the decays, as well as in the combined processes. A 
detailed numerical analysis has been performed. The dependences of the asymmetries 
on the C2HDM parameters are studied taking into account the theoretical constraints, 
experimental lower bound on the charged Higgs mass from B — > X s j and the constraint 
on the p parameter. 

The calculations have been performed with the help of the packages FeynArts and 
FormCalc. For that purpose, a new model file for FeynArts has been created and corre- 
sponding fortran drivers for FormCalc have been written. 

In the allowed parameter space of the C2HDM parameters, which is severely con- 
strained by vacuum stability, perturbativity, unitarity, lower charged Higgs mass bound 
and Ap, the studied CP violating asymmetries cannot be greater than ~3 %. This is in 
contrast with a similar study performed within the MSSM [13]. There the CPV asymme- 
tries can be larger by an order of magnitude. However, after having taken into account 
the relevant branching ratios and cross sections, the measurability of the studied CPV 
asymmetries in the C2HDM and in the MSSM at the LHC has roughly the same statistical 
significance, which is maybe not big enough for a clear observation at the LHC. At the 
SLHC with a design luminosity bigger by a factor of ~ 10, such a measurement would be 
worth of being performed. 

A The Ap constraint 

The p parameter is defined by the ratio of the neutral and charged currents at vanishing 
momentum transfer. In the C2HDM, at tree level p = my(c^m|) = 1 which is in 
perfect agreement with the experimentally measured value. This relation may be spoiled 
by radiative corrections if they are large. 
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For the extra contributions Ap to the p parameter from the additional scalars Higgses 
in the C2HDM, we have implemented in our code the expression given in [30] in terms of 
physical masses and elements Rjk of the rotation matrix (2.8). The analytic expression for 
Ap can be split into two contributions: 
i) Higgs-Higgs contribution (HH): 2 



^w(o) - cos 2 e w Aff (o) = ^ E [< t sin - cos P R irf + R h} F *p( M H± > M j) 



^[(sin/3% - cos 0R j2 )R k 3 ~ (sin/Jiifc! - cos/3^ 2 )i? j3 ] 2 F Ap (Mf, M fc 2 )] , (A.l) 
k>j 



with 



F Ap (mlml) = l(ml + ml) 



,2^,2 m 2 

lo e-2' 



mf — 777,2 



(A.2) 



ii) Higgs-ghost contribution (HG): 



a 2 r 

4w(0)-cos 2 %Aff (0) = -|_^g[cos/3% +sin/3^ 2 ] 2 x 

j 

(3Fa p (m1, M 2 ) - 3F Ap (M^, M 2 )) +3F Ap (Mf v , M 2 ) - 3F A/3 (M|, M 2 )] . (A.3) 



From this contribution one has to substract the SM Higgs contribution of a mass Mq [30] . 
The choice of Mq is usually consistent with the fit analysis. As a default value we take 
M = 120 GeV. 



B Expressions for A^ 



In order to work with the parameter set (2.13), we need the explicit relations between the 
parameters of the Higgs potential and the physical Higgs masses and rotation angles. It is 
straightforward to derive the Aj, 7 = 1,2, 3, 4, 5 as functions of our input parameters (2.13) 



2 The relevant couplings are given explicitly in appendix B of [30]. 



[29]: 



^ r .2 J2 t\ ,r2 i / „ „ „ i _ „ \ 2 ji 2 i / „ „ „ „ „ "\2 ji 7 r2 „2 , ,2] 



C /3 



^[ciC 2 M H o + (cis 2 s 3 + sic 3 ) M H o + (cis 2 c 3 - sis 3 ) M H o - s^n 

1 



Ai 

A2 = -3— ^[sfc^M^o + (C1C3 - S^SzYMjjQ + (cis 3 + SlS 2 C 3 ) 2 MfjO - C^yU 



As = —^{ci8i[<%M*o + (s 2 s 2 - 4)M 2 H o + (s|c§ - ^)ML] 
+s 2 c 3 s 3 (cl - sl){M 2 H0 -M 2 H0 )} + \[2M 2 H+ - fi 2 } , 

3 2 1) 

A 4 = ^ [alM^j + c 2 2 s 2 3 M 2 H o + c 2 2 c 2 3 M 2 H o + /i 2 - 2M 2 H+ ] , 
Re ( A s) = ^[- s i M H? " 44M 2 H o - c 2 2 c 2 3 M 2 H o + fJ 2 } , 
Im(A 5 ) = l —^{cj3[cic 2 s 2 M 2 - c 2 s 3 (cis 2 s 3 + sic 3 )M 2 

+c 2 c 3 (sis 3 - cis 2 c 3 )M^o] + s l 3[sic 2 S2Mf I o + c 2 s 3 (cic 3 - sis 2 s 3 )M^o 
-c 2 c 3 (cis 3 + s 1 s 2 c 3 )M^ I o}} , (B.l) 

with sp = sin/3, eg = cos/3 and // is given by eq. (4.2). 

The expressions (B.l) are implemented in our numerical code. Note that there are 
limits of non-CPV in the C2HDM [29], corresponding to the following values of the a- 
parameters: 

Q 2 = ±7r/2, 

a 3 = dbr/2, 

a 2 = a 3 = 0. (B.2) 

At these limits, the elements -Ri 3 , R23 and R 33 of the mixing matrix (2.8) vanish and 
the expression (2.15) becomes unstable. As this case is a subject of a different model 
convention, namely it is already the CP conserving 2HDM, our numerical code would 
produce an error message and interrupt the evaluation. 



C Interaction Lagrangian 

For our calculations we have used the following part of the C2HDM Lagrangian: 

C.0.1 Interactions of two quarks with a gluon 

C qqg = -9sT%G™q k ^ qi , M = l,2,3, a = l,...,8, q = t,b, (C.l) 
where T a /2 are the Gell-Mann matrices and g s is the SU(3) strong coupling constant. 
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CO. 2 Yukawa interactions of the neutral and the charged Higgses 

£itH° = t(hf.P L + h*P R )tH° j = 1, 2, 3, 



j 



h t,j = — - icpRjsjht, h t 



y/2 \/2mwsp 
£>bbH° = KhbjPL + hfjP R )bHl j = 1, 2, 3, 



3 



J ' — — \71 ~~ \71 ' 

Awr± = *(yt-P-L + ybP R )bH + + 6(y fe P L + y t P R )tH- , 

y t = htcp, y b = hbS/3 , (C.2) 

with S3 = sin/3 and eg = cos/3; Pl = (1 — 7s)/2 and Pr = (1 + 7s)/2 denote the left 
and right projection operators, respectively, Pjfc with j, k = 1,2,3 are the elements of the 
mixing matrix 1Z, see eq. (2.8). 

CO. 3 Two quarks-W and two quarks-ghost interactions 

CtbW± = -^^ PlW » + h tl P L tW-); 

Ak?± = KiftPL + y b P R )bG + + b(y b P L + y t P R )tG~, 

~ t _ 9 m t ~ b = gm b ^ c ^ 

CO. 4 Triple scalar interactions with neutral and charged Higgses 

£-H°H+H- = —(fHOH+H-)jHjH + H~, 3 = 1> 2 >3, 

3 9 

(fH°H+H-)j = c ts[s}{h - A 4 - Re(A 5 )) + dp\ 3 ]Rji + 

sp[<%(\2 - A 4 - Re(A 5 )) + s}\ 3 ]R j2 + spcpIm(\ 5 )R j3 . (C.4) 

CO. 5 Neutral Higgs-charged Higgs-W and neutral Higgs-charged Higgs-ghost 
interactions 

£h°h+w- = - c p R j2 + iR j3 )[H°d»H+W»- - H*d»H-W» + \, 



j = 1,2,3, 



C H o H+G - = —(fHOH+G-^HjE+G- + h.c, j = 1,2,3, 
9 



(fH°H+G-)j = s p[ s }(^ + Re(A 5 )) + c^(2Ai - A 3 - A345) - iIm(A 5 )]P j i + 
c/3[-c}{\4 + Re(A 5 )) - s}{2\ 2 - A 3 - A345) - iIm(A 5 )]Pj 2 + 
[i(A 4 - Re(A 5 )) + (c| - 4)Im(A 5 )]P i3 , 
A345 = A 3 + A 4 + Re(A 5 ), j = 1,2,3. (C.5) 
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CO. 6 Interactions of a neutral Higgs with W-W, W-ghost and ghost-ghost 



D Implementation of the complex 2HDM into FeynArts and FormCalc 

At present, there are several software packages for deriving Feynman rules and doing cal- 
culations of particle processes on the market e.g. [16, 65, 67]. For the purpose of our work 
it is convenient to use the FA package. Up to now the FA package does not include a 
model file for the 2HDM with CPV. In the following sections we describe how we have 
implemented the new model file for calculations in the C2HDM into FA, and how we have 
extended the corresponding FC fortran drivers. 

D.l FeynArts model file for the C2HDM 

At the moment, the diagram generator FA recognizes three generic particle physics models: 
the SM, the MSSM and the 2HDM [16]. The information about the physics properties of 
each model (fields, their propagators and their couplings) is collected in the correspond- 
ing model file. The model files for these three models exist in two different varieties: 
based only on the electroweak subset or including in addition quantum chromodynamics: 
SM(QCD) .mod, MSSM(QCD) .mod, THDM(QCD) .mod. All couplings are expressed in terms of 
the parameters of the relevant Lagrangian. 

The 2HDM model implemented in FA is based on the Higgs potential (2.1), but with 
m i2 = A6 = A7 = and all other parameters are real. The latter constraints refer to a 
model with the particle content of the SM, but a physical Higgs sector analogous to the 
one of the MSSM, see e.g. [68]. Usually, this is the most commonly used version studied 
in literature, also called the CP-conserving 2HDM due to absence of complex parameters. 
(Note that there exist another public software package for calculations in the CP-conserving 
2HDM [65].) In order to generalize the existing 2HDM model file for the general complex 
case, which is described by the Higgs potential (2.1), one should recalculate the couplings 
in terms of the full set of parameters in eq. (2.1). 3 The new FA model file CTHDM.mod 
consists of 270 couplings, which includes all Higgs interactions: 

• Higgs- vector boson interactions extracted from the covariant derivatives (2.16) 

• Triple and quartic Higgs self-interactions from the scalar potential (2.1) 

• Higgs interactions with Fadeev-Popov ghosts 

3 In this case the fields and their propagators are not affected. 




gmwicpRji + spR j2 )H^W+W^ , j = 1, 2, 3 , 

-j(cpRji + spR j2 )[H°d»G + W»- - H$& t G-W + ], j = 1,2,3 

— MC3A1 + s /3 X 345 )R jl + 

8p{a 2 p \ 2 + c}X 345 )R j2 - spcpIm(\ 5 )R j3 ]H?G+G-, j = 1,2,3. 



(C.6) 
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• Yukawa interactions 

Some of them are obtained explicitly e.g. in [29, 69]. The rest of the couplings needed to 
complete the model file are not concerned and we have taken them from the already existing 
model file THDM . mod. These are: vector boson self-interactions, fermion- vector interactions, 
and with QCD couplings: gluon self- interactions, gluon-ghost and gluon-quark interactions. 
The model file can be used independently for diagram generation in the general complex 
case of the 2HDM and it is not necessarily related to further calculations. There is a rule 
included for a switch to the case Xq = A 7 = 0, which describes the 2HDM with a softly 
broken Z2 symmetry (2.10). In order to go back to the CP-conserving case one must also 
set rn\ 2 to zero and consider A5 as a real parameter. The complete CTHDM .mod is too lengthy 
to be printed explicitly in this note. It can be found and downloaded from the FA website: 
www .feynarts. de . 

D.2 FormCalc drivers for the C2HDM 

After the diagram generation with the new FA model file CTHDM. mod FC calculates the 
squared matrix elements with the help of Form [71] and the resulting expressions are 
translated into Fortran for the further numerical evaluation. For consistency, the Fortran 
drivers necessary for the initialization of the model have to be extended to include the new 
set of parameters of the Higgs potential, the relations between them and the constraints 
on them. We would like to note that in spite of the fact that the new FA model file is 
written for the most general complex case with scalar potential (2.1), the extension of 
the FC fortran drivers is made only for the case of a softly broken Z2 symmetry of the 
2HDM Lagrangian, described by the scalar potential (2.10). The existing initialization file 
modeljthdm.F is replaced by the new file model.cthdm.F, which defines all parameters of 
the potential (2.10), the theoretical and experimental constraints on them. Also the set 
of input parameters for the numerical evaluation in the C2HDM is defined therein. All 
implemented expressions are listed and discussed in the present paper. 

Since the triple and the quartic scalar couplings have complicated and lengthly ex- 
pressions, during the calculation of a given process they are replaced by generic couplings 
named cS(i, j ,k) and qS(i, j ,k,l). Furthermore, all triple and quartic scalar couplings 
are evaluated in the model.cthdm.F file. 

D.2.1 Inputs 

Based on the case with a softly broken Z2 symmetry of the 2HDM Lagrangian described 
by the scalar potential (2.10) we work with the following set of input parameters: 



• 


MhO 


= M H o 


- mass 




• 


MHH 


= M H o 


- mass 


of # 2 ° 


• 


MHp 


= M H+ 


- mass 


of the charged Higgs 



4 Note that here we assume one Higgs doublet to couple with only up-type fermions and the other one- 
only with down-type fermions or the so called 2HDM type II. 
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• TB = i>i/i>2 - the ratio of the VEVs 

• rml2 = Re (777.12) ~~ real part of the 77742 parameter of the Higgs potential 

• alpl = ai - mixing angle 

• alp2 = «2 _ mixing angle 

• alp3 = - mixing angle 

• rmO = Mo = 120 GeV: reference point to substract the SM Higgs contribution from 
A/7, see section A 

Once the input parameters are initiated in run.F, in model.cthdm.F the mass of the heav- 
iest neutral Higgs M H o is calculated by using eq. (2.15). Then M H o has to satisfy the 
following conditions: 

i) M^j > 0, which means that H3 is not a tachyonic mode; 

ii) M H o < M H o < M H o, see section 2.1. 

When an allowed value for M H o is obtained, the code proceeds with the evaluation of the 
Aj as given in appendix B. Finally, the triple and the quartic scalar couplings are calculated 
in model_cthdm . F . 

D.2.2 Constraints 

The file model_cthdm. F contains the following constraints, see section 2.5 and section 4.1: 

• Vacuum stability 

• Perturbativity and unitarity 

• Ap constraint 

All constraints can be switched on and off. In addition, the upper and the lower Ap bounds 
can be modified easily. 
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